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Abstract. - We investigate recurrence phenomena in coupled two degrees of freedom systems. 
| It is shown that an initial well localized wave packet displays recurrences even in the presence of 

coupling in these systems. We discuss the interdependence of these time scales namely, classical 
period and quantum revival time, and explain significance of initial conditions. 
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Characteristics of quantum systems, which exhibit chaos in classical domain, have posed 
interesting questions for researchers. In one of these systems, namely, hydrogen atom in mi- 
crowave field, the discovery of phenomenon of quantum dynamical localization [1] proved as 
a land mark in the young field of quantum chaos. Fishman et al. connected dynamical local- 
ization with Anderson localization [2] , and proved it to be a generic property of periodically 
driven quantum systems. Later, the phenomenon was experimentally observed [3]. In this 
paper, we develop analytical treatment for quantum recurrence phenomena [4, 5] in systems 
which may exhibit chaos in classical domain and thus establish the phenomena generic to 
these systems. 

Quantum recurrences originate from the simultaneous excitation of discrete quantum lev- 
els [6]. The existence of recurrences has been investigated in atomic [7-15] and molecu- 
lar [16-19] wave packet evolution. Study of some of the periodically driven quantum sys- 
tems [20-22], and two-degree-of- freedom systems such as stadium billiard [23] indicates the 
presence of quantum recurrences in higher dimensional systems as well. Recently it is proved 
that the quantum recurrences are generic to periodically driven system which may display 
chaos in their classical counterpart [24]. In this paper, we study the phenomena of dynami- 
cal recurrences in general higher dimensional system, and provide mathematical foundations 
to the phenomena by calculating recurrence times, using semiclassical secular perturbation 
theory. At these time scales, namely classical period and quantum recurrence time, an initial 
excitation produced in the system displays its full or partial reappearance . 

A quantum wave packet in its early evolution follows classical mechanics and it reappears 
after a classical period following classical trajectory. Later, following wave mechanics it spreads 
and collapses, however, the discreteness of the quantum world leads to its restructuring. We 
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show that, (i) the phenomena of dynamical recurrences occur in higher dimensional quantum 
systems provided that at least two of its degrees of freedom are coupled, (ii) Furthermore 
we explain that, (a) the nonlinearity of the uncoupled systems, and, (b) the initial conditions 
on the excitation contribute to the classical and the quantum recurrence times occuring in 
the coupled multi-dimensional systems. (Hi) We also study the interdependence of these time 
scales for different kinds of dynamical systems. 

We write the general Hamiltonian of a system with any of its two degrees of freedom 
coupled [25,26], as 

H = H o (X) + XH c (X,0) (1) 

where, H is the Hamiltonian of the system in the absence of coupling, expressed in the action 
coordinates X = (Ii,I 2 ). Moreover, H c is the coupling Hamiltonian which describes coupling 
between X and is periodic in angle, 9 — (0\ 1 9 2 ), for nonlinear resonances in the system. The 
parameter A describes the strength of the coupling. We express the coupling Hamiltonian as, 

H c = J2H n (X)e me , (2) 

n 

where n — (ni,n 2 ). Whenever the frequencies f2 — 8H /dX obey the relation, n.fl = niOi + 
n 2 £l 2 = 0, resonances occur in the system. 

We may find classical chaos in the two degrees of freedom system, as a function of the 
coupling strength, whenever the two degrees of freedom are coupled. We consider that the 
coupling exists between I\ and I 2 . Within the region of resonance we find slow variations 
in action, hence, following the method of secular perturbation theory [27], we average over 
faster frequency, and get the averaged Hamiltonian for the Nth resonance, as H = Ho(I) + 
XV cos(Nip). Here, ip = 9\ — (M/N)02, I = h is the action corresponding to the angle 9\, V 
is the Fourier amplitude, and M and N are relatively prime integers [25,26]. The action I 2 
becomes the constant of motion. Moreover, Hq(I) expresses uncoupled averaged Hamiltonian. 
In case the coupling exists between Af degrees of freedom, we may apply the semi-classical 
secular perturbation technique (TV — 1) times to study the effect of most dominant degree of 
freedom. 

The energy of the excitation changes slowly when we produce it in the vicinity of TVth 
resonance of the higher dimensional system and it is narrowly peaked. Therefore, we expand 
the unperturbed energy, H (I), by means of Taylor's expansion around mean action I , and 
keep only the terms up to second order [28]. As a result, we express the Hamiltonian of the 
system, governing the evolution of the excitation in the vicinity of 7Vth resonance, as 

H = H (I ) + w(J - Jo) + |(J - hf + XV cos(Ncp). (3) 

Here, Hq(Iq) is the energy of the uncoupled system at the action / = 1$, and u> expresses 
the first order derivative of Ho with respect to action calculated at I and defines classical 
frequency of the excitation. The parameter ( defines the nonlinear dependence of the energy 
of the system on the quantum number. 

We introduce the transformation Nip = 29 and quantize the dynamics around the Nth 
resonance by quantizing the action [29,30], that is, 

h d _Nh d 

As a result, the Hamiltonian of the resonant system in the quantum domain reads 

N 2 Ch 2 d 2 h d - 
H - —l- W + + H (I ) + XV cos 29 . (5) 
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Hence, the quantum mechanical system evolves according to the Schrodinger equation, Htpk = 
Skipk- Here, is the wavefunction of the system in the region of a resonance and is, there- 
fore, required to fulfill periodicity condition, i.e. ipk{& + tt) = i/)(6). The corresponding 
eigenvalue £k defines the eigen energy of the system. With the help of transformation, 
ipk = 4>k exp (—2iu>8/(NC t h)), we map the Schrodinger equation on the Mathieu equation, 



if 
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+ a v — 2q cos 29 
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The 7r-periodic solutions of Eq. © correspond to even functions of the Mathieu equation 
whose corresponding eigenvalues are real [31]. These solutions are defined by Floquet states, 
i.e. <j) v {&) = e P v (0), where P u (0) is the even order Mathieu function. 

In order to obtain a 7r-periodic solution in ^-coordinate we require the coefficient of ip in 
the exponential factor to be equivalent to an even integer number, k. This condition provides 
the value for the index v as 



2A: 
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N N(h 

Therefore, we may express the eigen-energy of the system as 

£ft = — g — a v{k) {q) - — +H {I a ) , 



(9) 



(10) 



where, a v {q) is the Mathieu characteristic parameter [31]. 

In order to check this result we study the case of zero coupling strength, that is A = 0. 
In this case, the value for Mathieu characteristic parameter becomes a v (q = 0) = v 2 . This 
reduces the quasi-energy, £f., to Eq. © in the absence of coupling term, that is, A = 0, and 
in addition leads to define k as (7 — Io)/h. 

The initial excitation produced at 7 = Iq observes various time scales at which it reappears 
completely or partially during its evolution. In order to find these time scales at which an 
excitation in the quantum mechanical coupled higher dimensional system recurs, we employ 
the eigenenergy of the system [33,34]. 

These time scales, are inversely proportional to the frequencies lo^\ where lo^ = 

£k I dl( % \ when calculated at 7 = 7q. The index i describes the order of differ- 
entiation of the quasi energy With the increasing values for i we have smaller frequencies 
indicating longer higher-order recurrence times. 

and defines classical period of the higher dimensional coupled 



The time scale, T 



(i) 



system and is inversely proportional to u> 



(2) 



W In the absence of coupling a/ 1 ) reduces to 

ce time in the 
Here, we have 



u). The time scale, — and defines quantum mechanical recurrence time in the 

higher dimensional coupled systems. It has inverse proportionality with a/ 2 ' 



(2\)- 1 hd 2 £ k /dI 2 \ I=Io . 
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On substituting the value for £k from Eq. i|l(J|) in the defination for and lo^ 2 \ we 
obtain the classical period as, 

= [1 - M^]T^\ (11) 
and the quantum recurrence time for the coupled system as, 

Tj[ Q) = [1-MP]T^ Q) . (12) 

Here, the time scales Tq(= 2-klo^ 1 ) defines classical period and Tq\= 2-n:{^ / hC > )~~ 1 ) defines 
ianti 

Mo are given as, 



(cl) 

quantum revival time in the absence of coupling. The time modification factors Mo and 



2 \ lo 2 J (1 - a* 2 ) 2 



' (13) 



and 



where, 



M(Q> = lf^V 3 + ^ (14) 

2 ^ 2 y (1 - M 2 )3 ^ 



/i = ^7' (15) 

Equations l|ll|) and (|12fl express the classical period and the quantum recurrence time in the 
presence of coupling between two degrees of freedom as a function of the coupling strength 
A, nonlinearity (, associated with the uncoupled system and other characteristic parameters 
of the system. Analysis of Eqs. 113(1 and 1(14(1 leads us to the conclusion that the first terms 
of the modification factors Mi cV> and M^ depend quadratically on the coupling strength A, 
and, on the nonlinearity £ present in the initial uncoupled system. Whereas they are inversely 
proportional to the fourth power of frequency u> in both the cases. The second terms are 
function of /i, defined in Eq. ((15(1 . Hence, for the coupling strength A = we find = 
andTf>=T W) . 

case a: In the absence of nonlinearity, i.e for £ = 0, the time modification factor for the 
classical period Mo° 1 ^ and for the quantum recurrence time M^ vanish, which is evident 
from Eqs. ((13(1 and 1(14(1 . Thus, for linear coupled systems the quantum recurrences take place 
at infinite time, i.e — Tfi®' = oo. Nevertheless, the system displays recurrences after 

the classical period, T^ cl ' = T^ cl ' = 2-k/uj. Hence, in the coupled linear higher dimensional 
systems only classical periods exist. 

case b: For fi < 1, the time modification factors M^ and Mo are related as, 

AfJ°> = -M4 cl) = 3a, (16) 

where 

(17) 



i fxvc 

Thus, the classical period and the quantum recurrence time r| Q) in the presence of 



coupling are related with the Tq ; and Tq of the uncoupled system as 

gyM) rp{Q) 



AT (cl) 4T (Q) 



(18) 
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In the absence of coupling we find T { x cl) = T (c ° and T ( X Q) = T (Q) which fulfills Eq. fT5|t . 

The case may be achieved for a weak nonlinearity in the dominant degree of freedom of the 
coupled system, that is C << !■ The quantum recurrence time in coupled system T x ^ and 
in uncoupled system Tq®\ depend inversely on nonlinearity in the system, and are therefore 
much larger than classical periods T^ cl) and T (c ° in this case. 

As it follows from Eq. Ultjfl . the quantum recurrence time T x ^ reduces by 3aT , whereas, 

the classical period T x increases by olTq . Hence, we may conclude that the quantum 

dynamical recurrence time T x ^ reduces much faster than the classical period T x cl ^ in presence 
of a small nonlinearity in the system. 

Since, the modification factor a displays direct proportionality with the square of the 
nonlinearity parameter, £ 2 , in the system. In the asymptotic limit, i.e. for Q approaching 
zero, we get the result = Tq and T x = Tq = oo as discussed in case a. 

case c: In the presence of a relatively larger value of the nonlinearity parameter and/or 
highly quantum mechanical systems, we may consider fi > 1. We have the time modification 
factors related as 

= = -(3, (19) 



where 



P=-J^) ■ (20) 
2 \N 2 (h 2 ' 



Thus, the classical period T x and the quantum recurrence time T x ^ in the presence of 



A 

coupling are related with the Tq c1 ^ and of the uncoupled system as 



T (d) T (Q) 

rp{0l) ~ rp(Q) 



(21) 



The time modification factors for the classical period and the quantum recurrence time 
vanish as (3 reduces to zero. Therefore, Eqs. (|llfl and (|12|) provide us the asymptotic result, 
that is, T ( x cl) = T (c/) and T ( X Q) = T (Q) and we find that Eq. (EH holds. 

The parameter /3, which determines the modification both in classical period and in quan- 
tum recurrence time, is inversely dependent on the fourth power of the Planck's constant K. 
Hence, for highly quantum mechanical cases we note that /3 approaches zero and the times of 
recurrence remain unchanged. 

The classical period is inversely proportional to the classical frequency, a/ 1 -*, which is 
controlled by the initial excitation energy. We note that: (i) In presence of no nonlinearity 
(case a), only classical period exists as discussed above. Hence as frequency approaches zero, 
the classical period becomes infinity which is the case of an open system; (ii) For small value 
of the nonlinearity parameter in case b, as frequency becomes very small the classical period 
T^ c ^ and quantum recurrence time T x [32] changes following a square law dependence on 
coupling strength, A, as we find from Eqs. (|16|l and (|1 T|) . However, for larger values of the 
frequency the times T x and T x remain unchanged; (Hi) For relatively large value of the 
nonlinearity parameter, as in case c, the situation is opposite. For smaller value of frequency 
we find larger value for N, therefore, the modification factor reduces to zero. Hence, from 
Eqs. (HHJ and QQJ, we find that t[ c1) = T (c/) and T [ X Q) = T (Q) . However, for larger values of 
the frequency, we note that T [ x cl) and T [ X Q) vary following square law behavior as a function 
of modulation strength, A. 
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Fig. 1 - The change in recurrence phenomena for the wave packet originating from the center of 
a resonance: (a) We plot the square of auto-correlation function, C 2 , as a function of time, t, for 
a Gaussian wave packet. The wave packet is initially propagated from the center of a resonance 
in atomic Fermi Accelerator [32] and its mean position and mean momentum are zo = 14.5 and 
po — 1.45, respectively. We find that the wave packet experiences recurrence after each classical 
period. Thick line corresponds to our numerically obtained result and dashed line indicates quantum 
recurrence occuring for an oscillator in harmonic motion. The classical period is calculated to be 4-zr 
when the coupling strength, A = 0.3. As the coupling strength A becomes zero the evolution of the 
wave packet completely changes, (b) We plot square of the auto-correlation function for the wave 
packet in the absence of coupling, i.e. A = 0. We find that now the wave packet experiences collapse 
after many classical periods. Later, it displays recurrence at quantum recurrence time Tg . The 
inset displays the short time evolution of the wave packet comprising many classical periods for A = 0. 



For an excitation produced at the center of a resonance the quantum evolution of the 
system is characterized by Eq. iJSJ under the consideration of 6 w 0. Hence, the secular 
Hamiltonian of the coupled system for exact resonance case is 

Equation IL'l't describes the Hamiltonian of a harmonic oscillator. Hence, this provides an 
evidence that if the excitation originates initially from the center of a resonance it will expe- 
rience recurrences after each classical period, as shown in Fig. 1(a). The quantum recurrence 
time now is — Tq = oo. Thus the evolution is the same in every dynamical system 

in the case of exact resonance. The oscillator frequency is N^/(XV. Thus in case the cou- 
pling strength, A, reduces to zero the harmonic oscillator behavior disappears. The system 
now possesses only the classical period Tq°^ and the quantum recurrence time of the 

uncoupled system, we find this behavior in Fig. 1(b). This effect provides us information 
about level spacing around the center of a resonance [35] as well. For harmonic oscillator the 
spacing between successive levels is equal, hence we conclude that for quasiencrgy levels of 
the Floquet operator belonging to the center of resonance the level spacing is always equal. 

We conclude that in higher dimensional systems a coupling results in modifying the recur- 
rence times available for the uncoupled systems so far as the dynamics is considered close to 
a nonlinear resonance. The phenomena is helpful to improve the efficiency of the Recurrence 
Tracking Microscope [36]. The current results may help to identify [5] quantum acceleration 
modes (QAM) [37] and understand the dynamics associated. Moreover the suggested treat- 
ment may reveal the understanding of the dynamics of Bose Einstein Condensates in coupled 
systems [38]. 
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